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To describe the metal insulator transition in the new oxide TlS^CoOs we investigate the elec- 
tronic structure of its high temperature tetragonal phase by local density band structure (LDA) and 
model Hartree-Fock calculations. Within LDA we find a homogeneous metallic and ferromagnetic 
ground state, but when including the strong Coulomb interaction in the 3d shell more explicitly 
within the Hartree-Fock approximation, we find an insulating state of lower energy that exhibits 
both spin and orbital order. The instability of the metallic state towards the insulating one is driven 
by orbital degeneracy and a near degeneracy in energy of the states of intermediate (s=l) and high 
(s=2) spin. We also interpret our results in terms of a simple model. 
PACS numbers: 71.30.+h, 71.20.-b, 75.50.-y 

Introduction 

Six years ago, the observation of colossal magneto resistance in doped manganites has lead to renewed interest in 
transition metal oxides, see [Q for a review. These materials often exhibit complex phases with coexisting magnetic 
and orbital Q order and with Jahn- Teller like distortions. The origin of this complexity is the subtle interplay between 
different types of local Coulomb interactions within their 3d electronic orbitals. 

The Co 3+ analogue TlSr 2 Co05 of non superconducting TlSr 2 Cu05 is a recent addition to this class of materials. It is 
of perovskite structure with a first order transition at « 310 K from a tetragonal and metallic phase with ferromagnetic 
Curie-like susceptibility at high temperature to an orthorhombic and insulating phase with two inequivalent kinds of 
Co ions at low temperature j3j. Below T w 150 K antiferromagnetic long range order sets in. Special interest in this 
new oxide is due to its quasi two dimensional (2D) character and, from the point of view of applications, it may be 
useful that the metal insulator transition in it occurs at room temperature. 

In this paper we explain the 310 K metal insulator transition of TlS^CoOs in terms of an instability of its high 
temperature electronic structure. We show that the homogeneous ground state can lower its energy by breaking 
orbital degeneracy and by acquiring simultaneously both orbital and spin order. 

Method of analysis 

(i) Band structure calculation 

We start with a density functional band structure calculation in the high-temperature phase using the local density 
approximation (LDA) for the exchange and correlation potential, see for a review. Above the transition temperature 
of 310 K, the crystal structure is tetragonal with P4/mmm space group (see Fig. 1). As input parameters of the 
calculation, we used experimentally determined lattice parameters of a = 0.375 nm, c = 0.877 nm and relative 
strontium and apex oxygen (01) positions (that are not fixed by the space group) of z = 0.2903 and z — 0.2330, 
respectively j|. We employed a recently developed full-potential nonorthogonal local-orbital (FPLO) minimum basis 
scheme [gj that imposes no shape restriction on the potential. The calculation was scalar-relativistic, with the spatial 
extent of the basis orbitals controlled by a confining potential of (r/ro) 4 that was optimized with respect to the total 
energy. The thallium {5s, 5p, 5d, 6s, 6p}, strontium {4s, Ap, 5s, 5p, 5c?}, cobalt {3s, 3p, 4s, Ap, 3d} and oxygen 
{2s, 2p, 3d} orbitals were treated as valence orbitals, while the lower lying orbitals were treated as core states. The 
inclusion of Tl {5s, 5p, 5d}, Sr {4s, Ap} and Co {3s, 3p} states in the valence states was necessary to account for 
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non-negligible core-core overlaps. The O 3d states were taken into account to increase the completeness of the basis 
set. The results of this LDA calculation will be discussed in detail further below. 



(ii) Tight-binding Hartree-Fock calculation 



To take into account the strong Coulomb interactions in the Co d shell more accurately we could have used the LDA 
+ U approach described in detail in [[|. Instead, and in order to gain a more direct insight into the mechanism of the 
metal insulator transition at work, we included the Coulomb interactions in a minimal tight binding Hamiltonian of all 
cobalt 3d and oxygen 2p orbitals in the C0O2 plane and used the results of the previously described LDA calculation to 
determine the hopping parameters of this model (see for similar work). According to Koster-Slater type symmetry 
considerations |^|, the overlap between d- and p-orbitals may be parametrized in terms of two parameters which we 
denote as t Bg and t t2g , respectively. Due to the perovskite structure, there is no direct Co-Co hopping and we assume 
absence of direct 0-0 hopping. With these approximations, the transfer matrix elements for the planar Co-0 bond 
in x-direction are as follows 
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The superscript x for the p-orbitals indicates that they belong to the oxygen at \ae x and the corresponding matrix 
for the Co-0 bond in y-direction is obtained by formally interchanging x <-> y in the above table. The results of 



the band structure calculation determine also the Co d shell crystal-held parameters e,, i = {xy,yz, zx, z 



y 2 } 



and the on site energy e p of the oxygen p orbitals (we neglect the splitting of the oxygen p-states due to tetragonal 
distortion). Tetragonal symmetry at the Co site leaves e yz = e zx degenerate. The sum of on-site energies and of 
hopping amplitudes in the x- and y-directions defines the noninteracting part of a tight-binding model 



H — H t b + Hcoulomb 

H tb = ]T t l0 {R-R')a\ s {R)a. hS {R') 



(2) 



where the orbital indices i, j range from 1 to 11 (1 ... 5 are the Co 3d orbitals) and s denotes the spin component. To 
complete the model, we add the local Coulomb energy at the Co sites as expressed in tight-binding orbitals 
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To express the Vjjjy coefficients we use the rotational symmetry of the Coulomb matrix Vij^i in terms of spherical 
harmonics Y km |f 



V, 



77117712^1712 



t—l 

^ 2k + 1 J 

k=0 J 



47T / 

dndrz-^R^Rjirz) / d^ m J 2n ,Y km / dQ 2 Y 2 * m2 Y 2n2 Y^ 



= 25 (_)mi+n 2 p(k) 

fe=0,2,4 



2 2 k 




2 2k 
-mi n\ m 



F V» =A + -C, = 49B + 7C, = —C 

5 35 



2 2 k 

-712 TU2 m 



(4) 



and transform V to V in standard cubic orbitals given by 



In Rd(r) is the radial part of the 3c? wave function, {A, B, C} are Racah parameters and we used formula C.16, p. 

] . Please note that expression (Q) is rotationally invariant and that all Coulomb interactions were included. 
Although we used the full expression for Vijki in our calculation, the Racah parameter B is rather small and when it 
is taken to be zero we find a much simpler expression (in the cubic basis) 

( v ijki) B=0 = C (SijSki + SuSjk) + (A + C) S ik Sji (6) 

Usually, one denotes the diagonal elements as U = Van = A + 3C and the Coulomb repulsion between different (i ^ j) 
orbitals Vijij = A + C is smaller. However, in our context, it is more convenient to express the Coulomb energy in 
terms of the total electronic charge N = J2i n i an d total electronic spin stot = J2i &i P er site: 

{H C ouio m b) B=0 =CY1 P i P i - C z2si-Sj + — 4 — nin 3 

ij ij ij 

^J H [D-(s tot f] +u^- (7) 
with rij = a\ A ai^ + at ^4 , P { = a*, 1-0*4 and D = ^ P}Pi , 



with u = A + C/2, Jh = C and where = indicates that terms that merely renormalizc the chemical potential were 
ignored. In the last line we retained only the diagonal contribution of the pair hopping terms. 

Treating the Coulomb interaction in an unrestricted HF approximation JD| that allows for all spin and charge 
conserving correlations gives the following rcnormalization of the on-site Hamilton matrix 

l,m &' l,m 

where ij denote the orbital indices of Co 3d and the density matrix 

1 occ 

k,a 

is calculated from the normalized wave functions ^^(k) with band index a and momentum k (Nk corresponds to the 
number of k points). The total energy is given by 

E =< H t b > + < H Coulomb > 



where 
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is the mean (unrenormalized) kinetic energy, and 
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is the interaction energy. 



Results 



(i) Full-potential local-orbital method 



We performed two band structure calculations, one spin symmetric (see Fig. 2) and the other one allowing for 
spin polarization (Fig. 3). Both solutions are metallic but the magnetic one is energetically preferred by 0.54 eV per 
formula unit. The density of states (DOS) shows a high degree of covalency such that part of the magnetic moment 
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(calculated with all the overlap contributions) sits on oxygen (mo — 0.2/ib) giving a total moment m = 2.1/is 
(mco = 1.9/x.b). The corresponding occupation numbers are nc = 7.2 and no = 5.1. 

The band structure of the nonmagnetic solution together with the DOS for an easier identification of the structures 
is shown in Fig. 4; the size of the dots included in the band structure in Fig. 4 symbolizes the relative cobalt 3d- weight 
in the band. Evidently, the five bands crossing the Fermi level have predominantly Co 3d character. They hybridize 
quite strongly with 2p bands of the in-plane 0(2) located at about 4 to 6 eV binding energy. The 2p bands of the 
other oxygens 0(1) and 0(3) are nonbonding and located in between Co-3d and 0(2)-2p. The Co-3<i bands have 
only a small dispersion in z-direction (with the exception of Co 3d z 2 ) confirming the 2D character of the compound 
under consideration. To determine the energetic order of the 3d orbitals we used the bands at the V point plus 
their predominant orbital character fl4|| and found, in increasing order: d xy , d zx / yz , d z 2, d x 2_ y 2. Due to tetragonal 
symmetry, there is exact degeneracy between the d zx and the d yz orbital in the LDA calculation. 

The density of states of the ferromagnetic solution (Fig. 3) indicates a finite value for the majority spin part. 
From this we conclude that TlS^CoOs is not a half- metal and that it may conduct electric current even in the spin 
disordered phase. An orbital analysis (that is not documented here) shows that for majority spin only 3d x 2_ y 2 is 
partly occupied, while all the remaining majority bands (3d xy , 3d zx / yz , 3d z 2) are completely below the Fermi level. 
For minority spin, 3d xyi 3d zx / yz and 3d z 2 are partly filled, whereas 3d x 2_ y 2 is nearly empty, ft is difficult to interpret 
such an itinerant ferromagnet in an ionic picture. But its magnetic moment m = 2sub > 2ub indicates a state 
between intermediate (s = 1) and high spin (s = 2). 

(ii) Tight-binding Hartree-Fock calculation 

We used the nonmagnetic LDA solution of Fig. 4 to fix the parameters of the Hartree-Fock model (the magnetic 
solution contains the same information, except that the spin up and spin down bands are shifted against each other) . 
By comparing the five relevant bands of the FPLO result (Fig. 4) with the nonmagnetic HF solution we estimated 
the crystal field parameters to be (in eV): e xy = —1.0, e zx / yz = —0.5, e z 2 = —0.2, and e x 2_ y 2 = 0.5. The bandwidth 
determines the transfer terms t eg =1.9 and t t2cj — 1.4. The Racah parameters B and C are fixed to their ionic values 
for trivalent Co as determined by infrared spectroscopy (B — 0.06, C = 0.46) [jL5|, whereas A = 2 was chosen as a 
typical value for Co. The choice of e p depends on A due to the mean field shift of the 3d level and we used e p = 11 
such that the position of the oxygen levels (visible as the lower part in Fig. 6) coincides with the 0(2)-2p position 
in LDA. In fact, the position of the oxygen bands does not change very much for different Hartree-Fock solutions. 
The tight-binding bands resulting from these parameters are shown in Fig. 5 and agree reasonably well with the 
nonmagnetic LDA band structure of Fig. 4. 

Concerning the values of the the crystal field parameters, we note that due to tetragonal symmetry, there is exact 
degeneracy between the d zx and d yz orbitals in the LDA calculation. However, the values e xy = —1.0, e zx / yz = —0.5 
seem to violate " standard lore" || according to which the elongation of the octahedra in the z direction should shift 
d X y to higher energies compared to d zx / yz . One should note, however, that the £i are approximately the energies at 
the r point whereas the "standard lore" is valid in the ionic picture and would correspond, in our case, to the center 
of gravity of the different bands. Due to the larger bandwidth of the "two dimensional" d xy band in comparison with 
the "one dimensional" d zx / yz bands the corresponding centers of gravity nearly coincide in Figs. 4 and 5. For the 
electronic structure of the metallic, ferromagnetic high temperature phase of TlSr2Co05 it is crucial that all three 
bands {d xy , d zx / yz } cross the Fermi level which is due to their widths being larger than their crystal field splittings, 
see Fig. 4. It remains an open question, however, how the difference in dimensionality and width among the {d xy , 
d zx / yz } bands affects the screening of the Coulomb interaction. This question was considered in detail pi in the 
context of Sr2Ru04, another layered perovskite. 

We now discuss the magnetic solutions of the HF approach. The homogeneous ferromagnetic one (see Table I) 
is metallic and 650 meV lower in energy than the nonmagnetic solution. The good agreement with the energy gain 
in LDA (540 meV) and similar occupation numbers in LDA and HF indicate that our parameter assignment is 
satisfactory. Allowing for different occupations of d zx and d yz in a chess board like pattern, we find a metastable state 
with orbital order that decays into a ground state with both orbital and spin order and which contains Co sites with 
two different configurations, ms = 3.11/^s (close to high spin) and rriA = 2.07/is (intermediate spin). The energy gain 
due to the combined orbital and spin order is 67 meV. The origin of the instability of the homogeneous ferromagnetic 
state is the orbital degeneracy of d zx and d yz and the near degeneracy of intermediate and high spin configurations. 
This instability occurs in a rather large parameter region near the values which were derived for TlS^CoOs. However, 
due to the large dimensionality of the parameter space of our HF model (^) , we did not perform a systematic study. 
In Fig. 6 we give the spectral density of the ferromagnetic and of the two chess board like ordered states. As we can 
see, any kind of order leads to a decrease of spectral density at the Fermi level, but only the orbital- and spin-ordered 
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solution is insulating. That is also visible in the HF band structure (Fig. 7). 

The chess board like superstructure found here does not correspond to the orthorhombic low-temperature phase seen 
experimentally with a 2:1 ratio of high spin and intermediate spin states ||. We also investigated this experimental 
superstructure of the C0O2 plane within the model Hartree-Fock approach but found no solution with lower energy 
than those with chess board order. The lattice degrees of freedom may have to be included into the model to obtain the 
correct pattern, because the nearest neighbor Co-0 distance for high spin is probably larger than the corresponding 
distance for intermediate spin (as can be concluded from the analogy to the famous Invar alloys FexNii-a: JlT]). The 
present investigation of the electron system alone can only indicate the instability of the high temperature phase but 
it is not able to predict the correct low-temperature crystal structure. 



Interpretation of results 



We now interpret the results of our LDA and HF calculations in terms of a simplified model. As a first step, we 
recall that there is only d <-> p hopping and neither direct d <-» d nor p <-> p hopping. This and the fact that only cobalt 
d states are at the Fermi level allows us to extract an effective d «-> d hopping by eliminating the oxygen orbitals in 
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where a e{x, y, —x, —y} denotes the direction of hopping and Ae is the offset between the Co and oxygen bands. 
Actually, due to the high degree of d *-* p hybridization the second order expression (^) is certainly not sufficient to 
provide correct numbers for t e JJ but nonetheless it should give the correct matrix structure. Just like the original 
hopping the effective t e JJ is also anisotropic and orbitally dependent, with d zx and d yz electrons hopping in the x 
respectively y direction and forming one dimensional bands. Using the Coulomb energy of Eq. ([?]) this provides us 
with a simplified model that involves only d orbitals: 
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where jjj and u are renormalized values of Jh and U due to the elimination of the oxygen orbitals. But the precise 
amount of renormalization is difficult to calculate and beyond the scope of the present discussion. 

The above five band model must be simplified further to extract the relevant degrees of freedom. We first note 
that the crystal field parameters derived in the previous section are such that in the ionic case (t e JJ — 0) the three 
configurations of Fig. 8 are lowest in energy. We further highlight the orbitals closest to the Fermi level and which 
we believe to be itinerant by bold arrows. This suggests a minimal model with three itinerant electron species, 
namely d} x2 _ y 2 and the two degenerate minority spin bands d\ x and d yz . This physical picture is also supported 
by the ferromagnetic HF solution where only d 2 ,,2 is partly occupied among all the majority spin bands and the 

x y 

occupation of d^ 2 and d^ 2 _ 2 in Table I is indeed small (the occupation of d xy , however, deviates quite strongly from 
unity). The low energy sector responsible of the metal insulator and spin transition should be the competition between 



s = 1 



xy U, zx Tj, yz |, 3z 



| or x <-» y 



(11) 



s = 2: xyU, zx^\, yz f, 3z 2 - r 2 f, x 2 - y 2 | 



In other words, xy |J, is only a spectator orbital while {zx t, yz t, 3z 2 — r 2 f } provide a total spin of 3/2 and the exact 
degeneracy between d zx and d yz is kept as one of the driving mechanisms of the transition. The spin transition is 
then due to the competition between a down spin electron zx J, or yz J, and an up spin electron x 2 — y 2 f. Combining 



the spins of {zx |, yz f , %z l 
as follows: 



T} into an effective spin 3/2 degree of freedom, we can represent this competition 



5 



{Slzxl} or {Slyzl} & {S\, x 2 - y 2 T} (12) 

2 2 2 

where we used Hund's coupling to exclude misaligned spins. To describe these qualitative ideas more precisely, we 
propose the following model: 

H = e-£al s (R)a 3<s (R)- £ a],.^ aU#) ( 13 ) 

R,s R,R',s i=1 

-2j H Sr (%, 3 ~ g R,2 ~ g R,l) + \ n R,i n li,3 +jJ2^R' 

R R,i,j R,R' 

where the orbitals are numbered as {1,2,3} «-> {zx,yz,x 2 — y 2 } and is an anisotropic and orbitally dependent 
hopping matrix. Model ( p"3[ ) is formulated in an extended phase space in comparison to ([lO]) since |53,zx ]) and 

\S\ , yz 1) do not exist in the 5 band model. But those unphysical states are at a high energy due to Hund's coupling 

2 

The origin of the antiferromagnetic exchange coupling j should be the virtual superexchange of the singly 
occupied states that were excluded from our model and the coupling (^s^ 3 — sg 2 — makes sure that the 

spins of electrons in the orbitals i — 1,2 are antiparallel to the spin 3/2 vector and vice versa for i = 3. A detailed 
estimate of the model parameters in ([l3]) is beyond the scope of the present paper, but it is clear that we have to 
consider the range of parameters u ^> ju ^> £, t,J. The condition jg 3> t suppresses virtual hopping processes for 
antiparallel nearest neighbor spins £>g, whereas they remain possible for parallel spins. In the limit of jn ^> s,t the 
fermions are spin polarized, double occupation of the same orbital is automatically excluded by Fermi statistics and 
we therefore dropped the local double occupation .D-term in equation ( |l3| ) . The above model still admits high energy 
processes (involving energy costs of u and jn) that must be integrated out to obtain a true low energy model. 



To argue for the minimal model (13) we compare some of the possible two site clusters in second order perturbation 
theory (Fig. 9) (for simplicity we use identical transfer amplitudes t and restrict ourselves to leading terms). We see 
that for a ferromagnetic spin arrangement there is a competition between an orbitally ordered state s = 1 (realized for 
e > t 2 /u) and a mixed orbital and spin ordered state (e < t 2 /u) which is also realized in the model HF approach. Next, 
we see from the mixed spin cluster in the figure that hopping processes are influenced by the relative spin orientation. 
In this cluster, there is a competition between ferromagnetic order favored by a gain in delocalization energy ~ 
and antiferromagnetic order favored by a gain in magnetic energy ~ j. The ferromagnetic delocalization energy is 
expected to dominate j and therefore neighboring high spin and intermediate spin states should order ferromagnetically 
(see |^]) and it is easy to see that neighboring high spin states (s = 2) have an antiferromagnetic exchange in our 
simplified model (|l3|). 

There is a similarity between the model we propose here and the Zener double exchange model, ]l9| , except that 
our model has flipped spins for some of the orbitals to describe the spin transition and we also have a total of 3 
species of electrons per site at jn 3> t. For sufficiently small values of u there will be, by analogy with the double 
exchange model, a ferromagnetic and metallic phase because charge transport is possible for parallel spin orientations, 
while for u — > 00 the system is insulating. So we expect a rich phase diagram of this reduced model as a function of 
its parameters {e, i, u, j} with phases of mixed magnetic and orbital or spin order and which may be metallic or 
insulating as a function of its parameters. 



Conclusions 



Using LDA band structure calculations for the high temperature tetragonal phase and a HF approach in a minimal 
model of the perovskite plane, we found states that are lower in energy than the homogeneous and ferromagnetic 
state. More specifically we found that the homogeneous state is unstable towards orbital order and spin state dis- 
proportionation. We propose this instability to be the driving mechanism of the metal to insulator transition. A 
dynamical spin disproportionation above the transition temperature would be compatible with Mossbauer data on 
57 Fe-doped TlS^CoOs that suggest the existence of two inequivalent magnetic sites f2(i|| . Based on our calculations, 
we also proposed a simplified model with only three states per site at jjj 3> t plus an extra spin of S = | and which 
we argue to have a rich phase diagram as a function of its model parameters. 
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Figure Captions 

Fig. 1: Crystal structure of the high temperature tetragonal phase. 
Fig. 2: DOS of the nonmagnetic solution. 
Fig. 3: DOS of the ferromagnetic state. 

Fig. 4: Nonmagneticc LDA band structure. The relative cobalt 3d weight of the bands is symbolized via black dots in 
the bandstructure. For comparison we also included the DOS (see Fig. 2) at the right hand side of the band structure. 
The broad line at the Fermi level between X and M is comprised of 3 bands (d xy , d zx , and d z i) which are nearly 
degenerate. 

Fig. 5: Band structure of the nonmagnetic HF solution. 

Fig. 6: Comparison of spectral weights for three magnetic HF solutions. 

Fig. 7: HF band structure of the orbital- and spin-ordered solution. 

Fig. 8: The lowest ionic configurations of Co 3+ in TlSr 2 Co05. Configurations a and b correspond to intermediate spin 
states and are degenerate (d zx ,d^ z ) whereas c is a high spin state (d? x2 _ y2 )- The bold arrows correspond to itinerant 
states, the remaining states are localized and are combined into an effectice spin 3/2 (middle row). The right row 
introduces a notation. 

Fig. 9: Perturbation theory for configurations of a two site cluster. 



Table I: Occupation numbers of HF solutions. 
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